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The r-mode frequencies of the Tolman VII solution for the slowly rotating non-barotropic approx- 
imation within the low frequency regime are estimated. The relativistic correction to Newtonian 
r-mode calculations is shown as function of the tenuity -H and is shown to be significant only for 
very compact neutron stars. 



I. INTRODUCTION 

. The possibility that the fluid motion of rotating neutron stars is unstable to the emission of gravitational radiation 
' has motivated numerous recent studies of the so-called "r-mode" instability Although this is a relativistic effect, 
, arising from perturbations of the spacetime, the majority of available work has been focused on Newtonian and post- 
Newtonian approximations. In the first-order (in the angular velocity Vl) relativistic treatment the eigenfrequencies 
are real and the modes are determined by one ordinary differential equation ||^ in the low frequency regime for 
^ non-barotropic |^ stars. This singular eigenvalue problem then gives rise to a continuous spectrum, but can also 
admit a discrete r-mode solution |^(At higher frequencies the r-modes exhibit complex eigenfrequencies, losing their 
singular structure.) It has been further shown ||] that if the mode occurs inside the star, then it is associated with a 
' diverging velocity perturbation. It has been shown that the r-modes of low-index realistic polytropic models exhibit 
this unphysical feature, while constant density models do not j^, suggesting that the r-mode instability might be a 
' curiosity of unphysical models. It is more likely, however, that it is the mathematical description of the problem that 
is failing, not the underlying physics. Taking this view, we present a straightforward way to estimate the relativistic 
correction to Newtonian r-mode frequencies for realistic equations of state. Calculations are then peformed in the 
slowly rotating counterpart to the Tolman VII solution - the exact solution of Einstein's fleld equations that is currently 
\^ , the best fit to a neutron star |^ . We show that the maximum relativistic correction to Newtonian r-mode calculations 
' based on this model is <^ 23% for entirely causal solutions. 

Q II. THE EIGENVALUE PROBLEM 

cr| 

\^ , The static spherically symmetric metric in "curvature" coordinates Ip[ is 

ds^ = ds\ -I- r'^dB^ -I- r^sin^{e)dcl)^, (1) 
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■ where dp is a static Lorentzian 2-surface with coordinates {t,r). The source of (|^) is considered to be a perfect 
] mathematical fluid (of isotropic pressure p{r) and energy density p{r)) generated by radial streamlines of constant r. 
To first order in the angular velocity (17), a rotating star can be described by the stationary axisymmetric metric 

ds^ = -e'^'^Ut^ + e^^'^^dr^ + r^dd^ + r"^ sir? (d) dcf^ - 2w{r)r^ sin^{9)dtd<P, (2) 
where the functions zy(r) and A(r) are identical to those in (|l|). The function w{r) satisfies Hartle's equation 0] 
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4 .Jo'-z:^) + 731^(0 = 0, (3) 



subject to the boundary condition 



1 = r! = [w{r) + -^]r=R, (4) 

where j = e ' ' 2 * ' and w{r) = Vl — w{r). The axial modes for non-barotropic stars are determined by Kojima's 
equation ||] 

(a - W(r))(e-^W/^1' -ArTTT{r)h'^ - (8^T(r) - 4^ + iii±ll)/i„) + 16a^T(r)/io + 0{n'') + ... = 0, (5) 
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where = 4-, 

dr ' 

Tir)=p{r)+p{r), (6) 

and ho is a perturbation variable {I and m are spherical harmonic indices). Note that we have included higher-oreder 
corrections to this description. The problem is singluar at when w{rs) — a, the eigenvalue, which is related to the 
mode frequency [a) by 

Newtonian r- modes have a — 1. The fluid velocity perturbation (ua) vanishes at the boundary defined by the vanishing 
of the isotropic pressure. It is related to the perturbation variable by 

2rnw(r) \{r}-u(r) 
- - ^ ' re ^ ho. (8) 



2mw{r) ~ l{l + 1)((T + mVL) 



Kojima's equation (||) admits a continuous spectrum from the centre (w(0)) to the boundary {w{R)) and may 
also admit discrete r-mode solutions p^ . The singular point corresponding to these solutions cannot be inside the 
star, as the velocity perturbation would diverge In this scenario, therefore, eigenvalue solutions are only physical 
if W{R) < a < 1, whereas many realistic equations of state yield a < w{R). It is likely that this unphysical behaviour 
arises from the neglect of higher order terms near the r-mode, since the r-mode instability appears even in Newtonian 
theory These higher-order terms can only be significant when the 16 a tt T(r)/io term in (5) is small. For realistic 
equations of state, p(r) monotonically decreases through the solution, reaching a minimum at the boundary, which 
is defined by p{R) — 0. Hence, the eigenvalues must correspond to a value of W{r) close to the boundary of the 
solution. Therefore a useful approximation for computing the relativistic correction to Newtonian r-mode calculations 
for realistic equations of state is that a ~ w{R) . 



III. THE TOLMAN VII SOLUTION 

The Tolman VII solution jl^ is one of only two known physically reasonable exact spherically symmetric perfect 
fiuid solutions of Einstein's equations where the energy density vanishes at the boundary . The associated metric 



dsi 



sin(ln( 



( Jl - 5 ^ + 3 ^ + ^ - 5/6 + ^2(^ -5^ + 3 ^r'dr' 
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where C is a number chosen to give zero energy density at the boundary and is given by 




12 73/2 _ 23 + 4 V3 v/7 (7 - 2) 



o4 arctan((;) 



where 



^ 373/2 - 6 V7 + tan(ln(2))^/3v/7 (7 - 2) 

^ -3 73/2 tan(ln(2)) + 6 ^7 tan(ln(2)) + ^3^7 (7 - 2) ' 

^ — j^, and 7 is the tenuity (total radius to mass ratio, -p-). The solution is causal ||l^ for 7 >r 
pressure is not a "elementary" function. In contrast, the energy density is simply 



(10) 



(11) 



3.707. The isotropic 



where pc 
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p(0-Pc(l-e'), (12) 
Despite its simplicity, this density profile (and the corresponding binding energy and moment of 



inertia determined from it) is in good agreement with realistic equations of state for neutron stars with M > I.2M0 
0, even in the star's intermediate regions. Another appealing feature of this solution is that it does not have a phase 
transition (in p) at the boundary. Since the Tolman VII solution can reproduce several features that are independant 
of the neutron star equation of state, we take the view that using it is the best way by which to gauge the relativistic 
correction to Newtonian r-mode calculations. 
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IV. RESULTS AND CONCLUSIONS 



Numerical integrations of (3) in the Tolman VII solution give the plot shown in Figure 1. At the causal limit 
(subluminal adiabatic sound speed throughout the interior) the eigenvalue is a « w{R) « 0.77 which gives an 
estimate of 23% for the correction to the Newtonian r-mode calculation. The size of the correction then decreases 
monotonically with increasing 7. We find, therefore, that the difference between the r-mode frequencies of realistic 
slowing rotating non-barotropic neutron stars and those of Newtonian models in the low-frequency approximation, is 
only significant for very compact neutron stars. 
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Physically reasonable values of T 

1. 



generalized solution the energy density at the boundary is given by p{R) — 

for neutron stars are T — 1 + t where e <~ 10~*. In this paper, without any change to our conclusions, we take T 
Also note that the Tolman VII equation of state closely resembles a n=l polytrope at the boundary (the index decreases 
monotonically towards the centre of the star) . 

By causal we mean a subluminal adiabatic sound speed (v^ = ^ < 1). 

M. Ishak, L. Chamandy N. Neary and K. Lake, Phys. Rev. D 64, 024005 (2001) ( ^r-qc/0007d73| ). 

Portions of this work were made possible by use of GRTensorll. This is a package which runs within Maple. It is en- 
tirely distinct from packages distributed with Maple and must be obtained in dependently. The GRTensorll softw are and 
documentation is distributed freely on the World- Wide- Web from the address tittp: //grtensor .phy .queensu. ca 



3 




0.84 

0.83- 

0.82- 

0.81- 

0.80- 

0.79 

0.78 




3.8 



— [— 
4 



4.2 



4.4 



4.e 



4.8 



FIG. 1. Numerical integration of (3) in the Tolman VII solution. The solution is entirely causal for 7 >~ 3.707. An 
enlarged view of the compact region is shown. At the causal limit (subluminal adiabatic sound speed throughout the interior) 
the eigenvalue is a « w{R) 0.77 which gives an estimate of 23% for the maximum correction to the Newtonian r-mode 
calculation. Note that the size of the correction then decreases monotonically with increasing 7. 
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